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The electromagnetic properties of the baryon decuplet are calculated in quenched QCD on a 
20^ X 40 lattice with a lattice spacing of 0.128 fm using the fat-link irrelevant clover (FLIC) fermion 
action with quark masses providing a pion mass as low as 300 MeV. Magnetic moments and charge 
radii are extracted from the electric and magnetic form factors for each individual quark sector. From 
these, the corresponding baryon properties are constructed. We present results for the higher order 
moments of the spin-3/2 baryons, including the electric quadrupole moment E2 and the magnetic 
octupole moment MS. The world's first determination of a non-zero MS form factor for the A 
baryon is presented. With these results we provide a conclusive analysis which shows that decuplet 
baryons are deformed. We compare the decuplet baryon results from a similar lattice calculation 
of the octet baryons. We establish that the environment sensitivity is far less pronounced in the 
case of the decuplet baryons compared to that in the octet baryons. A surprising result is that 
the charge radii of the decuplet baryons are generally smaller than that of the octet baryons. The 
magnetic moment of the reveals a turn over in the low quark mass region, making it smaller 
than the proton magnetic moment. These results are consistent with the expectations of quenched 
chiral perturbation theory. A similar turn over is also noticed in the magnetic moment of the 
E*", but not for H* where only kaon loops can appear in quenched QCD. The electric quadrupole 
moment of the fl~ baryon is positive when the negative charge factor is included, and is equal to 
0.86 ±0.12 X 10~^ fm^, indicating an oblate shape. 

PACS numbers: 12.39.Fe, 12.38.Gc, 13.40.Em, 14.20.Dh, 14.20.Jn 



I. INTRODUCTION 

The study of the electromagnetic properties of baryons 
provides valuable insight into the non-perturbative struc- 
ture of QCD (see Refs. P, i, i, i, [| for recent reviews). 
Baryon charge radii and magnetic moments provide an 
excellent opportunity to observe the non-analytic be- 
havior predicted by chiral effective field theory (xEFT). 
Since these are inherently non-perturbative properties of 
hadrons, first-principles calculations on the lattice are es- 
sential for our understanding of hadronic structure, and 
indeed there has been much progress in this direction, 
mainly for the nucleon and pseudoscalar states (see @ 
for a review). For decuplet baryons, however, there has 
been very little progress since Ref. Q which appeared al- 
most 15 years ago. However, renewed interest appeared 
recently ^8, .9]. 

The Adelaide group has been investigating the electro- 
magnetic structure of hadrons for several years now. In 
Refs. [13, mi) we presented a novel method for determin- 
ing the strange quark contribution to the nucleon's elec- 
tromagnetic form factors, the results of which were later 
confirmed by an improved analysis of old experimental 
data [l^ and new data from parity violating experiments 
at JLab [l^. This was followed by an in-depth study in 
quenched QCD of the electromagnetic properties of the 
octet baryons [l3| . Of particular interest was an observed 
environmental isospin dependence of the strange quark 
distributions in A" and More recently, we performed 
an investigation into the pseudoscalar and vector meson 



electromagnetic form factors [l5|. Here we determined 
that the p"*" meson has a negative quadrupole moment, 
indicating that the p meson is oblate. 

In this paper we continue our study of the electromag- 
netic structure of hadrons and present a quenched lattice 
QCD calculation of the electromagnetic form factors of 
5f7(3)flavor decuplet baryons. From these form factors 
we determine magnetic moments, charge and magnetic 
radii and present results for the electric quadrupole and 
magnetic octupole moments. 

On the lattice, decuplet baryons are stable as a result 
of the unphysical large quark masses that are used in 
present calculations and the finite volume of the lattice. 
Decay to a pion and an octet baryon is forbidden by en- 
ergy conservation. However, stability of decuplet baryons 
is common to most hadronic models. In this sense, lattice 
results provide a useful forum in which the strengths and 
weaknesses of various models may be identified. The lat- 
tice results also provide access to observables not readily 
available with present experiments such as the higher- 
order multipole moments of the Q~ which is stable to 
strong interactions. 

An examination of decuplet baryon structure in lat- 
tice QCD enables one to study new aspects of non- 
perturbative quark-gluon dynamics. In analyzing the re- 
sults we make comparisons within the baryon decuplet 
and with the octet results which provide insights 
into the spin dependence of quark interactions. 

The E2 and M3 moments accessible in spin-3/2 sys- 
tems provide insights into the shape of the decuplet 
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baryon ground state. These higher-order moments also 
have the potential to discriminate between various model 
descriptions of hadronic phenomena. 

To put our results into perspective, we compare 
our calculations with experimental measurements where 
available, and with the predictions of Quenched Chiral 
Perturbation Theory (QxPT). 

The decuplet baryon interpolating fields used in the 
correlation functions are discussed in Sec. Ill Al The ex- 
traction of baryon mass and electromagnetic form factors 
proceeds through a calculation of two and three-point 
correlation functions. These are discussed in Sec. IIIBI 
The two and three point functions for decuplet baryons 
are discussed in Sees. Ill CI and IIIDI Throughout this 
analysis we employ the lattice techniques introduced in 
[Tj], and these are summarized in Sec. IIIII In Sec. IIVI we 
outline the methods used in our analysis of the lattice 
two and three point functions. Our results are presented 
and discussed in Sec. IVl and summarized in Sec. IVII 



II. THEORETICAL FORMALISM 

A. Interpolating Fields 

The commonly used interpolating field for exciting the 
A++ resonance from the QCD vacuum takes the long 
established [Iq, ^3 form of 



for the A"*" composed of two u quarks and one d quark 
has the form 



A+ / N 



1 



^abc 



V3 



2 (u^°(a;)C7^d^(2;)) u%x) 

+ {u^''ix)Cjf,u\x)) d^ix) 1 .(2.2) 



Other decuplet baryon interpolating fields are obtained 
with the appropriate substitutions of u{x), d{x) — > 
d{x) or s{x). The interpolating field for S**^ is 
given by the symmetric generalization 



3 *^ 



{u^''{x)Cj^d'>{x)) s%x) 
{s'^^{x)Cjf,u\x)) d%x) 1 .(2.3) 



The S'C/(2)-isospin symmetry relationship for E* form 
factors 



-1*-)- 



(2.4) 

may be easily seen in the E* interpolating field by noting 
e"''^ {s^''{x)Cjf,u''{x))d'{x) = 



^abc f^.Ta 



u^''{x)C-i^,s^{x))d%x). (2.5) 



A++/ X 



abc 



u^"(x)C7^u^(x)) u''{x). 



(2.1) 



Unless otherwise noted, wc follow the notation of Sakurai 
[l^ . The Dirac gamma matrices are Hermitian and sat- 
isfy {7^,7,y} = 2(5^^, with cr^^. = [7^,71^]- C = 7472 is 
the charge conjugation matrix, a, 6, c are color indices, 
u{x) is a u-quark field, and the superscript T denotes 
transpose. The generalization of this interpolating field 



B. Correlation functions 

Two-point correlation functions at the quark level are 
obtained through the standard procedure of contracting 
pairs of quark fields. Considering the A+ correlation 
function at the quark level and performing all possible 
quark field contractions gives the two-point function as 



T 



+2ST' tr 



tr 



257 tr 



7^ CS'^'''' C 7^ sx 



(2.6) 



where the quark-propagator S'^J" = T yu'^{x)u"' (0)j and 

similarly for other quark fiavors. 5'C/(3)fiavor symmetry 
is clearly displayed in this equation. 

In determining the three point function, one encoun- 
ters two topologically different ways of performing the 
current insertion. Figure [1] displays skeleton diagrams 
for these two insertions. These diagrams may be dressed 
with an arbitrary number of gluons. Diagram (a) illus- 



trates the connected insertion of the current to one of 
the valence quarks of the baryon. Diagram (b) accounts 
for the alternative time ordering where the current first 
produces a disconnected qq pair which in turn interacts 
with the valence quarks of the baryon via gluons. 

The number of terms in the three-point function is four 
times that in Eq. (|2.6p . The correlation function relevant 
for a A+ three-point function is 
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FIG. 1: Diagrams illustrating the two topologically diflterent 
insertions of the current within the framework of lattice QCD. 



where 



_ ^abc a'b'c 

3 



ASf lu CS™' C 7p Sf + ASf 7. CSj""' C -f^ 



AST 



ASf 7. CSj''''C^^ Sf + iSf 1. CS^'^'C^^ Sf + iSf 



j.CS^'Cj^S'T 
C Sj'''' C 7^ S^ 
7.C5™'C7^^r 



+2ST tr 
+2Sf tr 
+2Sf tr [• 



7^ CS'^^'' C S^" 

j,csr'''cj,,ST 



2ST tr 



2ST tr 



2ST tr 



luCSj''^ Cjf, s^ 
'^^CS^^'Cj^^S^ 



If C S'^'''' C 7p S^ 



^^CS™ C^fj^S^ 
'^,CS™'C-f^ST 

+ 2Sr tr '^.CSl'^'C-i^ST 



2ST tr 



-11 rl a 1 ^ 



q—u, d, s 



4Sr i.CS^'Cj^Sr +4Sr i^CSrC-ff^S^- +4S^, 



■Tbb' 



'j^CS^^'Cj^vS:^ 



+2ST tr 7, CS^'^'Cj,, ST + 2ST tr 7, CSrd^. + 2ST tr 7, CSf" C^^ S^ 



Tbb' , 



yTbb', 



ST {X2,x,,0) = V S^^X2,x^) 7a ^^^^ (xi, 0) , 



(2.7) 
(2.8) 



denotes the connected insertion of the probing current to 
a quark of charge e,. Here we have explicitly selected the 
electromagnetic current. However, the present discussion 
may be generalized to any quark-field-based current op- 
erator bilinear in the quark fields. 

The latter term of Eq. (|2.7p accounts for the discon- 
nected quark loop contribution depicted in Fig. [TJ). The 
sum over the quarks running around the loop has been re- 
stricted to the flavors relevant to the ground state baryon 
octet and decuplet. In the S'?7(3)flavor hmit the sum 
vanishes for the electromagnetic current. However, the 
heavier strange quark mass allows for a nontrivial result. 
Due to the technical difficulties of numerically estimating 
for the squared lattice volume of diagonal spatial 
indices, these contributions have been omitted from lat- 



tice QCD calculations of electromagnetic structure in the 
spirit of QxPT, and we will also do so here. For other 
observables such as the scalar density or forward matrix 
elements of the axial vector current relevant to the spin 
of the baryon, the "charges" running around the loop do 
not sum to zero. In this case the second term of Eq. (|2.7p 
can be just as significant as the connected term 19, 20|. 

An examination of Eq. (|2.7p reveals complete symme- 
try among the quark flavors in the correlation function. 
For example, wherever a d quark appears in the corre- 
lator, a u quark also appears in the same position in 
another term. An interesting consequence of this is that 
the connected insertion of the electromagnetic current for 
A° vanishes. All electromagnetic properties of the A° 
have their origin strictly in the disconnected loop contri- 
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bution. Physically, what this means is that the valence 
wave function for each of the quarks in the A resonances 
are identical under charge symmetry. 

C. Two-point Green functions 

In this and the following subsection discussing correla- 
tion functions at the hadronic level, the Dirac representa- 
tion of the 7-matrices is used to facilitate calculations of 
the 7-matrix algebra. It is then a simple task to account 
for the differences in 7-matrix and metric definitions in 
reporting the final results using Sakurai's notation. 

The extraction of baryon mass and electromagnetic 
form factors proceeds through the calculation of the en- 
semble average (denoted ( • • • )) of two and three-point 
Green functions. The two-point function for decuplet 
baryons is defined as 

^g-.p..-p/3a ^ ^ \T{xZ{x)x^M) I ^ ) -(2.9) 

X 

I 



Here Zb{p) denotes the overlap associated with our 
smeared source. Zb (p) is associated with the sink which 
need not have the same smearing. 

Equations (|2.14p through (|2.17p provide four correla- 



Here ft represents the QCD vacuum, F is a 4 x 4 
matrix in Dirac space and a, (3 are Dirac indices. 
The subscripts a, r are the Lorentz indices of the 
spin-3/2 interpolating fields. At the hadronic level 
one proceeds by inserting a complete set of states 
I B,p,s ) and defining 



( n I xAo) I b,p,s) = Zb{p)J — u,{p,s), (2.10) 

V p 



where Zb represents the coupling strength of Xct(O) to 
the baryon B. Our use of smeared interpolators makes 
this momentum dependent. Momentum is denoted by 
p, spin by s, and Ua{p, s) is a spin- vector in the Rarita- 
Schwinger formalism. Ep = \JpP' + Af ^ and Dirac indices 
have been suppressed. Using the Rarita-Schwinger spin 
sum. 



(2.11) 



(2.12) 



(2.13) 



(2.14) 
(2.15) 
(2.16) 
(2.17) 



tion functions from which a baryon mass may be ex- 
tracted. All baryon masses extracted from the different 
selections of Lorentz indices agree within statistical un- 
certainties. The combination providing the smallest sta- 



^Wa(p, s)Ur(p,s) = 



7 • p -I- Af f 1 2p^p^ p„^r - Prla 



9<yr - 7;lalT 



2M r 3' ' 3Af2 

Kr, 



our usual definitions for F, 



^^■^2^0 J ' 2\0 

and p = (p, 0, 0), the large Euclidean time limit of the two point function takes the form 

{ G^f{t;p,U) ) = ZB{p)ZBip)—e-''^' tr [ F4 A,, ] , 



where 



(G-(.;.^,r.)) ^ ZBip)z-Bip)l^^[^)e-^^^, 

{Gf,^{t;p,U)) = ZB{p)Z-B{p)\§,(^^^^^y-^^\ 

{Gg^{t;p,r,)) = ZB{p)z^{p)\{^^^^^yr'^^\ 



{ Gf3^(t;p, F4) ) = ZBipmp) ^'^'^ ^'^^ ^ 



3 V 2i?p 



e 



r 
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tistical fluctuations is ( G^^{t\p,Ti) + G^^{t]p,Vi) ) 
and these results are presented in section [V] 

It should be noted that the spin-3/2 interpolating 
field also has overlap with spin-1/2 baryons. For the A 
baryons and fi^ this poses no problem as these baryons 
are the lowest lying baryons in the mass spectrum having 
the appropriate isospin and strangeness quantum num- 
bers. However, E* and S* correlation functions may have 
lower lying octet spin-1/2 components allowed by flavor- 
symmetry breaking, m„ = rud ^ m^. Therefore it is 
desirable to use the spin-3/2 projection operator [2l[ 



1 



(2.18) 

However, our precision results for baryon two-point func- 
tions give no indication of a low- lying spin-1/2 compo- 
nent being excited by the spin-3/2 interpolating fields, 
and conclude such excitations are negligible. 



in terms of the covariant vertex functions and in terms 
of the current matrix elements. The electromagnetic cur- 
rent matrix element for spin-3/2 particles may be written 
as 



(p',,s' I j^(0) |p,s) = 



(2.19) 

Here p and p' (s and s') denote the momentum (spin) 
of the initial and final states, respectively, and Ua{p,s) 
is a Rarita-Schwinger spin-vector. The following Lorentz 
covariant form for the tensor 



g^P { air 



02 



2M 



2Mb 



} , (2.20) 



D. Three-point functions and multipole form 
factors 

Here we begin with a brief overview of the results of 
Ref. [22], where the multipole form factors are defined 



where P = p' + p, q = p' — p and Mb is the mass of the 
baryon, satisfies the standard requirements of invariance 
under time reversal, parity, G-parity and gauge transfor- 
mations. The parameters ai, 02, ci and C2 are indepen- 
dent covariant vertex functions. 



The multipole form factors are defined in terms of the covariant vertex functions through the following Lorentz 
invariant expressions |22l] . 

GEoiq^) = (1 + {a, + (1 + T)a2} - ^t{1 + r) {ci + (1 + r)c2} , (2.21) 

^£2(g') = {ai + + T)a2} - i(l + r) {ci + (1 + r)c2} , (2.22) 

GMiiq') - (1 + Jr)ai-^T(l + r)ci, (2.23) 

GMsiq^) = ai-^(l + r)ci, (2.24) 

with T — — g^/(2Ms)^ (> 0). The multipole form factors Geo, Ge2, Gmi and Gm3 are referred to as charge (EO), 
electric-quadrupole {E2), magnetic-dipole (Ml) and magnetic-octupole (M3) multipole form factors, respectively. 

In a manner similar to that for the two-point function, the three-point Green function for the electromagnetic 
current is defined as 

( G^f^{t2,h;p\p;T) ) = ^ g-^p'.x-.g+<p'-p)..lp/3a ^ ^ | y (^a(^2)jnxi)x?(0)) (2.25) 



Once again, the subscripts a, r are the Lorentz indices of the spin-3/2 interpolating fields. For large Euclidean time 
separations ^2 — ti » 1 and ti » 1 the three-point function at the hadronic level takes the limit 

( G^f^ih,h;p',p;r) )^ ^e-^p'(*-*i)e-^''*^r'3" ( | | p',-s' )( p',s' \ j^^ \ p, s ){ p, s \x^, \ n ), (2.26) 

where the matrix element of the electromagnetic current is defined in (j2.19p . and the matrix elements of the interpo- 
lating fields are defined by Eq. (|2.10p . 

The time dependence of the three-point function may be eliminated through the use of the two-point functions. 
Maintaining the lattice Ward identity, which guarantees the lattice electric form factor reproduces the total charge of 
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the baryon at (f' = 0, provides an indispensable guide to the optimum ratio of Green functions. The preferred ratio 
of two- and three-point Green functions is 7] 



R„^r{p'J\T), (2.28) 



where we have defined the reduced ratio R^^ t-{p' iP] T)- Note that there is no imphed sum over a and r in Eq. (|2.27p . 
Also, the square root in Eq. (|2.27p spoils the covariant/contraviant nature of R^^r and no meaning should be attached 
to the location of the indices. We still prefer this notation due to the close connection with G^^''^ . 

Using our standard definitions for F given in Eq. (|2.12D and the Rarita-Schwinger spin sum of Eq. (|2.1ip . the 
multipole form factors may be isolated and extracted from the following combinations of Ra'^ r{p' ,P', T) 

QEo{q^) = i(i^l^(gl,0;^4)-f i?2S(<fi,0;r4)-fi?3''3('ri,0;r4)), (2.29) 



gE2{q^) = 2^^^^^±^(i?iS((fi,0;r4)+i?2'2(gl,0;r4)-2i?3'3(gi,0;r4)) 



(2.30) 



GMiiq^) = -l^rt^ {Ri\{qu0;r2) + R2^2iqi,0;r2) + R3\iqi,0;T2)) , (2.31) 

5 \qi\ 

GMsiq^) = -4 ^'^^^^^3^^^' (fll^(gl,0;^2) + fi2^2(gi,0;r2) - ^i?3^3(gi,0;r2)j , (2.32) 

where qi = (5,0,0). We note that smaller statistical uncertainties may be obtained for Qe2 by using the symmetry 

i?2S(gi,0;r4) -i?3S(gi,0;r4), (2.33) 

4 

in Eq. (|2.30p . Hence, we define an average R^^^g as 

Rtgiqi,0;T4) = i p?2'2(<fi,0;r4) +:R3'3(9i,0;r4)] . (2.34) 



With this definition the expression for QE2{q^) used in our simulations is 

GE2{q^) = 2 ^'^(^+/-^) (i?iS(qi,0;r4) - ELg(gi,0;r4)) . (2.35) 



III. LATTICE TECHNIQUES 

The three-point functions discussed in section |TT] are 
constructed usingthe sequential source technique out- 
lined in Refs. [tI. l23l. [23| . Our quenched gauge fields are 
generated with the C(a^) mean- field improved Liischer- 
Weisz plaquette plus rectangle gauge action using 
the plaquette measure for the mean link. The simula- 
tions are performed on a 20"^ x 40 lattice with a lattice 
spacing of 0.128 fm as determined by the Sommer scale 
j26| tq = 0.50 fm. This large volume lattice ensures a 
good density of low-lying momenta which are key to giv- 
ing rise to chiral non-analytic behavior in the observables 
simulated on the lattice [l^, llH • 

We perform a high-statistics analysis using a large 
sample of 400 configurations for our lightest eight quark 
masses. We also consider a subset of 200 configurations 



for our three heaviest quark masses to explore the ap- 
proach to the heavy-quark regime. A small sub-ensemble 
bias correction is applied multiplicatively to the heavy 
quark results, by matching the central values of the 200 
configuration sub-ensemble and 400 configuration ensem- 
ble averages at k = 0.12780. All tables display the raw, 
unbiased data for the first four kappa values. The first 
row of the n = 0.12780 results gives the results from the 
200 configuration sub-ensembles, the second gives the 400 
configuration ensemble results. The scaled results from 
the 200 configuration sub-ensembles are shown in the fig- 
ures. 

We use the fat-link irrelevant clover (FLIC) Dirac op- 
erator [27I which provides 0{a) improvement [2^. The 
improved chiral properties of FLIC fermions allow effi- 
cient access to the light quark-mass regime [1^, making 
them ideal for dynamical fermion simulations now under- 
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way [30|. For the vector current, we an C'(a)-improved 
FLIC conserved vector current We use a smeared 

source at to = 8, and a current insertion centered at 
ti = 14. Complete details are described in Ref. 
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Table [J provides the kappa values used in our simula- 
tions, together with the calculated n and decuplet baryon 
masses. While we refer to to infer the quark masses, 
we note that the critical value where the pion mass van- 
ishes is Kcr = 0.13135. 

We select k — 0.12885 to represent the strange quark 
in this simulation. At this k the ss pseudo scalar mass is 
0.697 GeV, which compares well with the experimental 
value of 2 m^ — m^ — (0.693 GeV)^, motivated by leading 
order chiral perturbation theory. 

The error analysis of the correlation function ratios is 
performed via a third-order, single-elimination jackknife, 
with the per degree of freedom (Xdof ) obtained via co- 
variance matrix fits. We perform a series of fits through 
the ratios after the current insertion at ti = 14. By ex- 
amining the Xdoi able to establish a valid window 
through which we may fit in order to extract our observ- 
ables. In all cases, we required a value of Xdof larger 
than 1.5. The values of the static quantities quoted in 
this paper on a per quark-sector basis correspond to val- 
ues for single quarks of unit charge. 

When extracting form factors from the lattice corre- 
lation functions via the ratios defined in Eqs. (|2.29p to 
p.32p in Sec. Ill D[ we employ the advanced analysis tech- 
niques outlined in detail in Ref. 14]. 

The following calculations are performed in the lab 
frame p = 0, p' = q — |q| x at \q\a = 27: /L^ with 
Lx = 20, the minimum nonzero momentum available on 
our lattice. While is dependent on the mass of the 
baryon, we find this mass dependence to be small. In- 
deed all form factors may be regarded as being calculated 



at 



= 0.230 ± 0.001 GeV^ where the error is 



dominated by the mass dependence of the target baryon. 
Where a spatial direction of the electromagnetic current 
is required, it is chosen to be the z-direction. 



IV. CORRELATION FUNCTION ANALYSIS 

A. Baryon Masses 

Figure [2] is a plot of the decuplet baryon masses along 
with the masses of the octet baryon from our previous 
calculation [Ti"!. We observe the 5'C/(3)fiavor limit at our 
sixth quark mass when = 0.485(3) GeV^. The mass 
of the ri~ is the mass of the A at the S'C/(3)fiavor limit, 
i.e., 1.732 ±0.012 GeV which differs from the experimen- 
tally measured value of 1.67 GeV by only about 3.6%. 
The higher value from the quenched simulation is in ac- 
cord with the expectations of quenched xEFT jsT], Is^] . 
The mass of the A baryon shows an upward chiral cur- 
vature as the becomes smaller. This behavior has 
already been discussed in Refs. [29l[3ll. [33j . 
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FIG. 2: Masses of the decuplet and the octet baryons at 
different quark masses. The octet baryon masses are the open 
symbols while the symbols of the decuplet baryons are filled. 
At the lowest quark mass the lowest point is the Nucleon, 
followed by A, E and the H. The decuplet baryon with lowest 
mass is the A followed by E* and H*. 



B. Form Factor Correlators 

The baryon form factors are calculated on a quark- 
sector by quark-sector basis with each sector normalized 
to the contribution of a single quark with unit charge. 
Hence to calculate the corresponding baryon property, 
each quark sector contribution should be multiplied by 
the appropriate charge and quark number. Under such 
a scheme for a generic form factor /, the A+ form fac- 
tor, /^+ , is obtained from the u- and d-quark sectors 
normalized for a single quark of unit charge via 

/4+ =2x ^ x/„ + l X (^-0 x/d. (4.1) 

Figure [3] depicts the electric form factor i?0 of the u 
quark in the A as a function of time at the SU (3)fiavor 
limit. The u and d quarks in the A states are identical 
as discussed in regard to Eq. (|2.7p . The straight lines 
indicate the fits which were selected using the x^^j con- 
siderations outlined in Ref. [l^. For light quark masses 
smaller than the strange quark mass, we fit the change 
in the form factor correlation functions from one quark 
mass to the next and add this to the previous result at 
the heavier quark mass. This provides significant can- 
cellation of correlated systematic errors and makes the 
selection of the fit regime transparent. 

Figure 2] shows the fitting of the electric form factor 
splitting for the A"*" between the eighth and ninth quark 
masses, the latter having m\ = 0.215(4) GeV^ The 
improvement of the plateau is apparent in Fig. 31 Still 
substantial Euclidean time evolution is required to obtain 
an acceptable x^of 
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TABLE L Hadron masses in units of GeV for various values of the hopping parameter, k. Pion masses are in GeV^ while the 
baryon masses are in GeV. 







A 


E* 




0.12630 


0.9960(56) 


1.999(9) 


1.908(10) 


1.815(11) 


0.12680 


0.8936(56) 


1.945(10) 


1.871(11) 


1.797(12) 


0.12730 


0.7920(55) 


1.890(10) 


1.834(11) 


1.779(13) 


0.12780 


0.6920(54) 


1.836(11) 


1.798(12) 


1.761(13) 


0.12780 


0.6910(35) 


1.845(10) 


1.807(11) 


1.770(11) 


0.12830 


0.5925(33) 


1.791(11) 


1.771(11) 


1.752(12) 


0.12885 


0.4854(31) 


1.732(12) 


1.732(12) 


1.732(12) 


0.12940 


0.3795(31) 


1.673(14) 


1.693(13) 


1.712(13) 


0.12990 


0.2839(33) 


1.622(16) 


1.659(15) 


1.695(13) 


0.13205 


0.2153(35) 


1.592(17) 


1.638(15) 


1.685(13) 


0.13060 


0.1384(43) 


1.565(18) 


1.620(16) 


1.676(14) 


0.13080 


0.0939(44) 


1.549(19) 


1.609(16) 


1.670(14) 


experiment 


0.0196 


1.232 


1.382 


1.531 




o 
cd 
6h 



o 



o 



16 18 20 22 24 
Euclidean Time 



28 




0.4 



16 18 
Euclidean Time 



FIG. 3; EG electric form factor of the u quark in the A at 
= 0.230(1) GeV'^ as a function of Euclidean time {t2) for 
ml = 0.4854 GeV^ the S'(7(3)flavor limit. The line indicates 
the fitting window and the best fit value. 



Tables HI] to IIVI list the electric form factors for all 
the decuplet baryons at the quark level for the eleven 
quark masses considered. In the tables, the selected time 
frame, the fit value and the associated x^^j are indicated. 
Table IVl provides collected results for the various decuplet 
baryons. 

The magnetic form factor Ml for the u quark sector in 
the A at the SU(3) limit is plotted in Fig.Oas a function 
of Euclidean time. Here the conversion from the natu- 
ral magneton, e/{2mB), where the mass of the baryon 
under investigation appears, to the nuclear magneton, 
e/(2mAr), where the physical nucleon mass appears, has 
been done by multiplying the lattice form factor results 
by the ratio tun / tub- In this way the form factors are 



FIG. 4: EO electric form factor of the u quark in the A at 
= 0.230(1) GeV^ as a function of Euclidean time {t2) 
at the ninth quark mass where ml = 0.215(4) GeV^ The 
correlator is obtained from the splitting between the ninth 
and eighth quark mass results. The line indicates the fitting 
window and the best fit value. 



presented in terms of a constant unit; i.e. the nuclear 
magneton. 

In Fig. [S] we present the Euclidean time dependence 
of the the Ml magnetic form factors of A calculated at 
the ninth quark mass where = 0.215(4) GeV^ using 
the splittings analysis. The onset of noise at this lighter 
quark mass is particularly apparent at time slice 20. 

Results for the quark-sector contributions to the Ml 
magnetic form factors of decuplet baryons are summa- 
rized in Tables |Vll to IVIIII While some of the Xdof 
some what large we note that neighboring regimes with 
acceptable Xdoi h&ve a variation in the central value that 
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TABLE II: Quark sector contributions to the electric form factor EO of the A at = 0.230(1) GeV'^. Sector contributions are 
for a single quark of unit charge. The fit windows are selected using the criteria outlined in Ref. ,14]. The quark contribution 
at the SU(3) limit when = 0.485(3) GeV^ provides the s quark contribution in Q~ . 



ml 










dA 




(GeV^) 


fit value 


fit window 


Xdof 


fit value 


fit window 


Xdof 




U.oUUt^O ) 


on 0/1 


1/17 


n &f\f\f 

U.oUUl^O ) 


on 0/1 

ZU — z4 


1/17 




n 7Q/1 

u. / y^t^D ) 


on 0/1 


u.yo 


n 7Q/1 f fi^ 


on 0/1 
zu — z4 


n QK 

U.yo 


u. / yzui^oo ) 


U. / oo\ t ) 


on 0/1 

ZU — z4 


u. / y 


n 7Sft/'7^ 


on 0/1 

ZU — z4 


n 70 
u. / y 




U. ( oo\o ) 


on 0/1 

ZU ~ z^ 


u.oy 


n 79i'^^9i^ 

U. ( 00(^0 J 


on 0/1 

ZU — Z4 


n pio 
u.oy 


U.DyiUl^oo ) 


\j. ( ( 6[p ) 


on 0/1 
zU — z4 


1.4D 


n 77Q('f^^ 


on 0/1 
zU — z4 


1 /ifi 


0.5925(33) 


0.764(7) 


20 - 24 


1.07 


0.764(7) 


20 - 24 


1.07 


0.4854(31) 


0.755(10) 


20 - 24 


0.53 


0.755(10) 


20 - 24 


0.53 


0.3795(31) 


0.744(11) 


17-20 


0.73 


0.744(11) 


17-20 


0.73 


0.2839(33) 


0.733(13) 


17-20 


0.79 


0.733(13) 


17-20 


0.79 


0.2153(35) 


0.725(16) 


17-19 


0.46 


0.725(16) 


17-19 


0.46 


0.1384(43) 


0.717(22) 


17-19 


0.23 


0.717(22) 


17-19 


0.23 


0.0939(44) 


0.693(33) 


17-19 


0.34 


0.693(33) 


17-19 


0.34 



TABLE III: Quark sector contributions to the electric form factor EO of S* at = 0.230(1) GeV'^. Sector contributions are 
for a single quark having unit charge. The fit windows are selected using the criteria outlined in Ref. [l^ . 



ml 




lis* or rfs* 










(GeV^) 


fit value 


fit window 


Xdof 


fit value 


fit window 


Xdof 


0.9960(56) 


0.804(6) 


20 - 24 


0.85 


0.759(9) 


20 - 24 


0.54 


0.8936(56) 


0.798(7) 


20 - 24 


0.65 


0.761(9) 


20 - 24 


0.50 


0.7920(55) 


0.792(8) 


20 - 24 


0.63 


0.763(10) 


20 - 24 


0.49 


0.6920(54) 


0.786(9) 


20 - 24 


0.56 


0.766(10) 


20 - 24 


0.48 


0.6910(35) 


0.774(6) 


20 - 24 


1.15 


0.752(8) 


20 - 24 


0.97 


0.5925(33) 


0.764(8) 


20 - 24 


0.91 


0.753(8) 


20 - 24 


0.79 


0.4854(31) 


0.755(10) 


20 - 24 


0.53 


0.755(10) 


20 - 24 


0.53 


0.3795(31) 


0.744(11) 


17-20 


0.89 


0.754(10) 


17-20 


0.22 


0.2839(33) 


0.733(12) 


17-20 


0.81 


0.754(11) 


17-20 


0.37 


0.2153(35) 


0.727(14) 


17-19 


0.15 


0.753(11) 


17-19 


0.21 


0.1384(43) 


0.719(18) 


17-19 


0.07 


0.753(12) 


17-19 


0.05 


0.0939(44) 


0.710(23) 


17-19 


0.22 


0.746(14) 


17-19 


0.63 



is small with respect to the statistical uncertainty. Re- 
sults for the various decuplet baryons are given in Ta- 
ble CXl 



This behavior has also been observed in recent lattice 
calculations [34|. Using this observation, together with 
the standard small expansion of the Fourier transform 
of a spherical charge distribution 



DISCUSSION OF RESULTS 



A. Charge radii 



It is well known that the experimental electric (and 
magnetic) form factor of the proton is described well by 
a dipole ansatz at small 



Ge{0) 



g2 > 0. 



(5.1) 



(5.2) 



we arrive at an expression which allows us to calcu- 
late the electric charge radius of a baryon using our 
two available values of the Sachs electric form factor 
(eis(Qii„): Qsm, namely 



i) 



12 



/ gg(0) 



(5.3) 
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TABLE IV: Quark sector contributions to the electric form factor EO of ^* at — 0.230(1) GeV^. Sector contributions are 
for a single quark having unit charge. The fit windows are selected using the criteria outlined in Ref. [13 ]. 













Us* 






fi t VP\ 1 1 1 P 


fit windnw 


Adof 


fi i" vf\ 1 1 1 p 


fi t wi Ti d nw 


y3 , 

Adof 


0.9960(56) 


0.765(10) 


20 - 24 


0.47 


0.809(8) 


20 — 24 


0.66 


0.8936(56) 


0.766(11) 


20 - 24 


0.47 


0.802(8) 


20 — 24 


0.58 


0.7920(55) 


0.767(11) 


20 - 24 


0.49 


0.795(9) 


20 — 24 


0.61 


0.6920(54) 


0.769(12) 


20 - 24 


0.50 


0.788(10) 


20 — 24 


0.58 


0.6910(35) 


0.753(9) 


20 - 24 


0.73 


0.775(7) 


20 - 24 


0.89 


0.5925(33) 


0.754(9) 


20 - 24 


0.66 


0.765(8) 


20 - 24 


0.75 


0.4854(31) 


0.755(10) 


20 - 24 


0.53 


0.755(10) 


20 - 24 


0.53 


0.3795(31) 


0.754(10) 


17-20 


0.19 


0.744(10) 


17-20 


1.19 


0.2839(33) 


0.754(10) 


17-20 


0.22 


0.734(11) 


17-20 


0.77 


0.2153(35) 


0.754(10) 


17-19 


0.03 


0.727(13) 


17-19 


0.09 


0.1384(43) 


0.755(11) 


17-19 


0.11 


0.720(14) 


17-19 


0.10 


0.0939(44) 


0.754(11) 


17-19 


0.66 


0.714(17) 


17-19 


0.19 



TABLE V: -EO form factor of the charged decuplet baryons for different values. The EO form factor of the A at the 



Sf/(3)flavor limit (to 


I = 0.485(3) GeV- 


^) provides the EO form factor of Q 








mi (GeV^) 


A++ 


A+ 


A- 








0.9972(55) 


1.601(10) 


0.803(5) 


-0.803(5) 


0.819(6) 


-0.789(7) 


-0.780(9) 


0.8936(56) 


1.589(11) 


0.794(6) 


-0.794(6) 


0.810(6) 


-0.786(7) 


-0.778(10) 


0.7920(55) 


1.577(13) 


0.788(7) 


-0.788(7) 


0.801(7) 


-0.782(8) 


-0.777(10) 


0.6920(54) 


1.566(16) 


0.783(8) 


-0.783(8) 


0.792(8) 


-0.779(9) 


-0.775(11) 


0.6910(35) 


1.545(11) 


0.773(6) 


-0.773(6) 


0.781(6) 


-0.766(7) 


-0.760(8) 


0.5925(33) 


1.527(14) 


0.764(7) 


-0.764(7) 


0.768(7) 


-0.761(8) 


-0.757(9) 


0.4854(31) 


1.509(19) 


0.755(10) 


-0.755(10) 


0.755(10) 


-0.755(10) 


-0.755(10) 


0.3795(31) 


1.487(22) 


0.744(11) 


-0.744(11) 


0.740(11) 


-0.747(10) 


-0.751(10) 


0.2839(33) 


1.465(26) 


0.733(13) 


-0.733(13) 


0.726(13) 


-0.741(12) 


-0.747(11) 


0.2153(35) 


1.451(31) 


0.725(16) 


-0.725(16) 


0.718(15) 


-0.736(13) 


-0.745(11) 


0.1384(43) 


1.433(44) 


0.717(22) 


-0.717(22) 


0.708(20) 


-0.730(15) 


-0.743(12) 


0.0939(44) 


1.386(65) 


0.693(33) 


-0.693(33) 


0.698(27) 


-0.722(19) 


-0.741(13) 



However to calculate the charge radii of the neutral 
baryons, the above equation cannot be used, due to the 
fact that in those cases Ge = 0- For the neutral baryons 
it becomes a simple matter to construct the charge radii 
by first calculating the charge radii for each quark sector. 
These quark sectors are then combined using the appro- 
priate charge and quark number factors as described in 
Sec. lIVIjj to obtain the total baryon charge radii. Indeed, 
all baryon charge radii including the charged states are 
calculated in this manner. 

Tables |X] to IXIII provide the electric charge radii of 
the decuplet baryons and the quark sector contributions. 
Figures [7] and [8] depict plots of the quark contributions 
to the decuplet charge radii. At the SU(3) limit (sixth 
quark mass) the quark contributions are identical in all 
cases as expected. 

In Figs, [ni and [TO] we compare the charge radii of the 



quark sectors in the decuplet baryons to those in the octet 
baryons at the ninth quark mass. From the figures it is 
evident that the contribution of the quarks is very much 
baryon dependent in the octet case, while for the decu- 
plet baryons there is much less environmental sensitivity 
to the individual quark contributions. More specifically, 
in the case of the u quark in the octet baryons, the charge 
radius decreases with the inclusion of the s quark, while 
such an influence of the s quark on u quark charge ra- 
dius is absent in the decuplet behavior. Furthermore, we 
note that the charge radius of the u quark distribution 
in the decuplet baryons is smaller than that in the octet 
baryons. 

From Fig. [9] it is evident that the charge radius of the u 
quark in the proton (up) is larger than that of the u quark 
in the A"*" (ua+)- In order to investigate this difference 
more accurately, we compute the ratios of charge distri- 
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TABLE VI; Quark sector contributions to the magnetic form factor Ml of A at = 0.230(1) GeV^. Sector contributions 
are for a single quark having unit charge. The fit windows are selected using the criteria outlined in Ref. [l3|. The quark 
contribution at the 5'f/(3)flavor limit when = 0.4854(31) GeV^ provides the s quark contribution in Q~ . 



ml (GeV^) 




y-A (mjv) 






dA (mjv) 






fit value 


fit window 


Xdof 


fit value 


fit window 


Xdof 




1 1 7'^^'9Pi^ 
i. 1 * oyzo f 


1 Q 9/1 






1 Q 9/1 
ly — z^ 


1 A ^ 




1 oni (9Q\ 


1 Q 94 

ry — 


1 1 A 


1 9ni ^'9Q^ 

1 .zun^zy ^ 


1 Q 94 

ry — z^ 


1 1 4 




1 9^ni'^^~> 


1 Q 94 


u.yo 




1 Q 94 

ry — Z4; 


u.yo 




1 9'\Pt^d^^ 


19 — 24 


U. ( o 




19 — 24 


n 7^1 


u.oy lui^oo } 


L.A^oyoZ f 


1 Q 9/1 


1 9f^ 


1 9/15i^''^9^ 


1 Q 94 

ly — z^ 


1 9Pi 
i .zo 


0.5925(33) 


1.269(41) 


19 - 24 


0.79 


1.269(41) 


19 - 24 


0.79 


0.4854(31) 


1.280(56) 


19 - 24 


0.31 


1.280(56) 


19 - 24 


0.31 


0.3795(31) 


1.301(64) 


17-21 


1.37 


1.301(64) 


17-21 


1.37 


0.2839(33) 


1.312(75) 


17-19 


1.14 


1.312(75) 


17-19 


1.14 


0.2153(35) 


1.309(89) 


17-19 


0.91 


1.309(89) 


17-19 


0.91 


0.1384(43) 


1.28(12) 


17-18 


1.26 


1.28(12) 


17-18 


1.26 


0.0939(44) 


1.11(21) 


17-18 


1.79 


1.11(21) 


17-18 


1.79 



TABLE VIL Quark sector contributions to the magnetic form factor All of E* at — 0.230(1) GeV^. Sector contributions 
are for a single quark having unit charge. The fit windows are selected using the criteria outlined in Ref. ,14]. 



ml (GeV^) 




ue* or cJe* (Miv) 






se* (Miv) 






fit value 


fit window 


Xdof 


fit value 


fit window 


Xdof 


0.9960(56) 


1.191(31) 


19 - 24 


1.07 


1.268(40) 


19 - 24 


0.51 


0.8936(56) 


1.216(35) 


19 - 24 


0.90 


1.275(42) 


19 - 24 


0.49 


0.7920(55) 


1.242(39) 


19 - 24 


0.82 


1.282(45) 


19 - 24 


0.47 


0.6920(54) 


1.264(45) 


19 - 24 


0.69 


1.289(50) 


19 - 24 


0.47 


0.6910(35) 


1.254(37) 


19 - 24 


0.94 


1.275(41) 


19 - 24 


0.69 


0.5925(33) 


1.272(45) 


19 - 24 


0.64 


1.278(47) 


19 - 24 


0.51 


0.4854(31) 


1.280(56) 


19 - 24 


0.31 


1.280(56) 


19 - 24 


0.31 


0.3795(31) 


1.297(62) 


17-21 


1.72 


1.289(60) 


17-21 


0.70 


0.2839(33) 


1.306(68) 


17-19 


1.88 


1.299(65) 


17-19 


0.04 


0.2153(35) 


1.305(76) 


17-19 


1.84 


1.309(69) 


17-19 


0.08 


0.1384(43) 


1.299(89) 


17-19 


3.44 


1.330(74) 


17-19 


0.44 


0.0939(44) 


1.25(12) 


17-18 


2.55 


1.303(88) 


17-18 


0.16 



butions of similar quarks in the octet to that in the de- 
cuple!. The uncertainty in the ratio (r^)(up)/(r^)(uA+) 
is calculated using the jack-knife method. Figures [11] 
and [13 depict the ratio of the quark contributions in the 
octet baryons to those in the decuplet baryons at the 
5?7(3)flavor limit and the ninth quark mass, respectively. 
In both cases, the doubly represented u quark contri- 
bution to the charge radius of octet baryons is larger 
than that in case of the singly represented octet quarks 
and decuplet quarks. At the SU(3) limit, all quarks take 
the strange quark mass, and hence one would expect the 
quark model picture to dominate. This suggests that 
the UA should have a broader distribution than the dis- 
tribution than that of Up due to hyperfine interactions. 



Our results contrast this prediction. The smaller charge 
radius of wa compared to that of Up also rules out any 
suggestion of a hyperfine attraction leading to ud diquark 
clustering in the nucleon or hyperon states [ssj . 

The charge radii of the various decuplet baryon states 
are shovifn in Figs. [131 HH and [15] as a function of m^. 
The charge radius of A~ is numerically equal to that of 
the with a negative sign. The charge radius of is 
taken as that of the A" in the S'f/(3)flavor limit, and is 
numerically equal to —0.307 ± 0.015 fm^. As our calcu- 
lations neglect the fl~ — > 'EP tt~ dressing, we anticipate 
our result to underestimate the magnitude. 

The decuplet baryon form factors are dominated by the 
net charge of the light quarks. For the A'^ the symmetry 
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TABLE VIII: Quark sector contributions to the magnetic form factor All of H* at = 0.230(1) GeV^. Sector contributions 
are for a single quark having unit charge. The fit windows are selected using the criteria outlined in Ref. T4]. 



mi (GeV^) 




sh* (mjv) 






or ds* (/iiv) 






fit value 


fit window 


Ador 


fit value 


fit window 


Aaoi 


0.9960(56) 


1.286(50) 


19 — 24 


0.45 


1.208(39) 


19 — 24 


0.84 


0.8936(56) 


1.289(52) 


19 — 24 


0.46 


1.231(42) 


19 — 24 


0.79 


0.7920(55) 


1.293(54) 


19 — 24 


0.46 


1.254(46) 


19 — 24 


0.76 


0.6920(54) 


1.297(56) 


19 — 24 


0.45 


1.273(51) 


19 — 24 


0.63 


0.6910(35) 


1.278(48) 


19 — 24 


0.48 


1.260(44) 


19 — 24 


0.68 


0.5925(33) 


1.280(51) 


19 - 24 


0.40 


1.274(49) 


19 - 24 


0.52 


0.4854(31) 


1.280(56) 


19 - 24 


0.31 


1.280(56) 


19 - 24 


0.31 


0.3795(31) 


1.285(58) 


17-21 


0.70 


1.293(60) 


17-21 


2.56 


0.2839(33) 


1.289(60) 


17-19 


0.04 


1.300(63) 


17-19 


2.66 


0.2153(35) 


1.293(62) 


17-19 


0.02 


1.303(66) 


17-19 


1.44 


0.1384(43) 


1.302(64) 


17-18 


0.73 


1.303(72) 


17-18 


2.36 


0.0939(44) 


1.301(67) 


17-18 


0.23 


1.313(81) 


17-18 


0.56 



TABLE IX: 


Ml form factor of the charged decuplet bary 


3ns for different 


ml values. 


The Ml form factor 


of the A at the 


Sf/(3)flavor 


imit {ml = 0.485(3) GeV^) 


provides the Ml form factor of Q. 








ml (GeV^) 


A++ 


A+ 


A- 


E*+ 






0.9972(55) 


2.35(5) 


1.17(2) 


-1.17(2) 


1.16(3) 


-1.22(3) 


-1.26(5) 


0.8936(56) 


2.40(6) 


1.20(3) 


-1.20(3) 


1.20(3) 


-1.24(4) 


-1.27(5) 


0.7920(55) 


2.46(7) 


1.23(3) 


-1.23(3) 


1.23(4) 


-1.26(4) 


-1.28(5) 


0.6920(54) 


2.51(8) 


1.26(4) 


-1.26(4) 


1.26(4) 


-1.27(4) 


-1.29(5) 


0.6910(35) 


2.50(6) 


1.25(3) 


-1.25(3) 


1.25(4) 


-1.26(4) 


-1.27(5) 


0.5925(33) 


2.54(8) 


1.27(4) 


-1.27(4) 


1.27(4) 


-1.27(5) 


-1.28(5) 


0.4854(31) 


2.56(11) 


1.28(6) 


-1.28(6) 


1.28(6) 


-1.28(6) 


-1.28(6) 


0.3795(31) 


2.60(13) 


1.30(6) 


-1.30(6) 


1.30(6) 


-1.29(6) 


-1.29(6) 


0.2839(33) 


2.62(15) 


1.31(7) 


-1.31(7) 


1.31(7) 


-1.30(7) 


-1.29(6) 


0.2153(35) 


2.62(18) 


1.31(9) 


-1.31(9) 


1.30(8) 


-1.31(7) 


-1.30(6) 


0.1384(43) 


2.56(24) 


1.28(12) 


-1.28(12) 


1.29(10) 


-1.31(8) 


-1.30(6) 


0.0939(44) 


2.22(43) 


1.11(22) 


-1.11(22) 


1.23(14) 


-1.27(10) 


-1.30(7) 



of the u and d quarks makes the form factors vanish. This 
is in contrast to the neutron where the three quarks are 
in mixed-symmetric states, giving rise to a non-zero form 
factor and charge radius. Charge radii of the neutral S* 
and S* are also close to zero and are dominated by the 
light quark sectors. 



The presence of the s quark as one moves from A to 
S* and S* reduces the magnitude of the charge radius as 
indicated in Figs. fT3l [Ml and [TSl By examining the ratio 
of the charge radii of the octet to decuplet baryons in 
Fig. [ini we observe that the octet baryons have a slightly 
larger charge radius than their decuplet counterparts. 



B. Magnetic moments 

The magnetic moment is provided by the value of the 
magnetic form factor at zero momentum transfer, = 
0, 



2Mf. 



(5.4) 



in units of the natural magneton, where Mb is the mass 
of the baryon. Since the magnetic form factors must be 
calculated at a finite value of momentum transfer, Q^, 
the magnetic moment must be inferred from our results, 
Qm{Q^), obtained at the minimum non- vanishing mo- 
mentum transfer available on our periodic lattice. We 
choose to scale our results from Gm{Q^) to Gm{0)- We 
do this by assuming that the dependence of the elec- 
tric and magnetic form factors are similar at the quark 
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TABLE X: Charge radii of the A baryons for different mi. Quark sector contributions for a single quark of unit charge are 
included. The charge radii of the A"'"+ are twice that of the A+ and the results for the A" are in QQCD. At the S'?7(3)flavor 
limit we find the charge radius of the Q,~ to be equal to —0.307(15) iia^. 



m^(GcV^) 


UA(fm^) 




rfA(fm^) 




A+(fm^) 


A"(fm^) 


0.9960(56) 


0.238(7) 




0.238(7) 




0.238(7) 


-0.238(7) 


0.8936(56) 


0.247(8) 




0.247(8) 




0.247(8) 


-0.247(8) 


0.7920(55) 


0.256(10) 




0.256(10) 




0.256(10) 


-0.256(10) 


0.6920(54) 


0.264(12) 




0.264(12) 




0.264(12) 


-0.264(12) 


0.o9iU(ooj 


O.Z(9(9j 








O.Z(9(9j 


— O.Z(9(9j 


0.5925(33) 


0.293(11) 




0.293(11) 




0.293(11) 


-0.293(11) 


0.4854(31) 


0.307(15) 




0.307(15) 




0.307(15) 


-0.307(15) 


0.3795(31) 


0.324(17) 




0.324(17) 




0.324(17) 


-0.324(17) 


0.2839(33) 


0.343(21) 




0.343(21) 




0.343(21) 


-0.343(21) 


{j.21oS[So) 


0.355(zd) 




0.355(26) 




0.355(26) 


-0.355(26) 


fl 1 QC/1 (A 

u. Loo^y^o ) 


U.o / U(^o / ) 




0.370(37) 




0.370(37) 


-0.370(37) 








0.410(57) 




0.410(57) 


-0.410(57) 


iAoL/hj Ai: Unarge 


radii of the S* baryons 


for different 


mi . Quark sector contributions for a single quark of unit charge are 


inplnHpH 

±X±^± LiV_l V_yV_X . 














m^(GcV^) 


Me* (fm^) 


5E.(fm2) 




E*+(frn2) 


S*"(fm2) 


E*"(fm2) 


0.9960(56) 


0.233(8) 


0.299(13) 




0.212(8) 


-0.022(3) 


-0.255(10) 


0.8936(56) 


0.242(10) 


0.296(14) 




0.224(9) 


-0.018(2) 


-0.260(11) 


0.7920(55) 


0.251(11) 


0.293(15) 




0.237(10) 


-0.014(2) 


-0.265(12) 


0.6920(54) 


0.260(13) 


0.289(16) 




0.250(12) 


-0.010(2) 


-0.270(14) 


0.69i0(ooj 


U.2(o[H)) 


V.Oll[lZ ) 




0.267(9) 


-0.011(1) 


-0.289(10) 


0.5925(33) 


0.292(12) 


0.309(13) 




0.286(11) 


-0.006(1) 




0.4854(31) 


0.307(15) 


0.307(15) 




0.307(15) 


0.000(1) 


-0.307(15) 


0.3795(31) 


0.324(17) 


0.308(16) 




0.330(17) 


0.006(1) 


-0.319(16) 


0.2839(33) 


0.341(20) 


0.309(17) 




0.352(21) 


0.011(1) 


-0.330(19) 




0.352(23) 


0.309(17) 




0.366(25) 


0.014(2) 


-0.338(21) 


0.1384(43) 


0.365(29) 


0.311(18) 




0.383(33) 


0.018(5) 


-0.347(25) 


0.0939(44) 


0.380(39) 


0.321(22) 




0.399(45) 


0.020(7) 


-0.360(32) 


TART.F, XTT- Hharcr 


e radii of the S* baryons for different values. Quark sector contributions for a sinj 


;le quark of unit 


charge are included. 














m^(GeV^) 


ss*(fm2) 




«a*(fm2) 




H*°(fm^) 


E*-(fm2) 


0.9960(56) 


0.291(16) 




0.227(11) 




-0.042(5) 


-0.269(14) 


0.8936(56) 


0.289(16) 




0.236(12) 




-0.035(5) 


-0.271(14) 


0.7920(55) 


0.287(17) 




0.246(13) 




-0.027(4) 


-0.273(15) 


0.6920(54) 


0.285(17) 




0.256(14) 




-0.019(3) 


-0.275(16) 




0.309(13) 




0.276(11) 




-0.022(3) 


-0.298(12) 


0.5925(33) 


0.308(14) 




0.291(13) 




-0.011(2) 


-0.302(14) 


0.4854(31) 


0.307(15) 




0.307(15) 




0.000(0) 


-0.307(15) 


0.3795(31) 


0.307(15) 




0.324(16) 




0.011(1) 


-0.313(16) 


0.2839(33) 


0.308(16) 




0.341(19) 




0.022(3) 


-0.319(17) 


0.2153(35) 


0.307(16) 




0.351(21) 




0.029(4) 


-0.322(17) 


0.1384(43) 


0.307(16) 




0.363(24) 




0.037(7) 


-0.326(19) 


0.0939(44) 


0.308(17) 




0.372(29) 




0.043(10) 


-0.330(20) 
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FIG. 5: Ml magnetic form factor of the u quark in A at 
= 0.230(1) GeV^ as a function of Euclidean time {t2) 
for = 0.485(3) GeV^ the 5'[/(3)-flavor hmit. The line 
indicates the fitting window and the best fit value. 
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FIG. 6: Ml magnetic form factor of the u quark in A at 
= 0.230(1) GeV^ at the ninth quark mass where m„ = 
0.215(4) GeV^. The correlator is obtained from the splitting 
between the ninth and eighth quark mass results. The line 
indicates the fitting window and the best fit value. 



FIG. 7: Charge radii of the u quark in A, E* and H* at 
different quark masses. The values for A are plotted at 
while that of the E* and H* are plotted at shifted for 
clarity. 
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FIG. 8: Charge radii of the s quark in E* and S* at different 
quark masses. The values for E* are plotted for shifted 
values for clarity. 



masses simulated herein. This is supported by experi- 
ment where the proton ratio ^if ^fLl-. — 1 for values of 
similar to that probed here. In this case 

5M(0) = |^g^,(0). (5.5) 

The strange and light sectors of hyperons will scale dif- 
ferently, and therefore we apply Eq. (|5.5p to the individ- 
ual quark sectors for all the decuplet baryons. Decuplet 
baryon properties are then reconstructed via 

QfAQ) = Qii{Q) + QlAQl (5.6) 



where I labels the light quarks and s labels the strange 
quark, and quark numbers and charges are included. 
Similar calculations are performed for the u and d sectors 
of the A. 

In Figs. [iTl [18] and [191 we display the quark sector 
contributions to the decuplet magnetic moments, which 
are equal in the SU (3)flavor limit (sixth quark mass) for 
single quarks of unit charge. Here we observe that the 
quark contributions in the S* and S* are very similar, 
which provides evidence that there is little environmental 
sensitivity. The chiral behavior of the ua is also very 
interesting. 
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FIG. 9: Charge radii of the u quark sector in the octet (open 
squares) and the decuplet (filled squares) baryons at the ninth 
quark mass where ml = 0.215(4) GeV^ 



FIG. 11: Ratio of charge radii of the quark sector contribu- 
tions in the octet/decuplet baryons at the S'(7(3)flavor limit 
where ml = 0.485(3) GeV^. 
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FIG. 10: Charge radii of the s quark sector in the octet (open 
squares) and the decuplet (filled squares) baryons at the ninth 
quark mass where = 0.215(4) GeV^. 



FIG. 12: Ratio of charge radii of the quark sector contribu- 
tions in the octet/decuplet baryons at the ninth quark mass 
where ml = 0.215(4) GeV^ 



C. Effective Moments 

In order to compare the decuplet moments with the 
octet moments from our previous lattice calculation , 
we construct so-called effective moments by equating the 
lattice quark sector contributions to the same sector of 
the SU(6)-magnetic moment formula derived from SU(6)- 
spin-flavor symmetry wave functions. 

The simple quark model formula from SU(6)-spin- 
flavor symmetry gives the magnetic moment of the pro- 
ton as 



_ 4 1 



(5.7) 



where fj.^ and fj,d are the intrinsic magnetic moments of 
the doubly represented u and singly represented d con- 



stituent quarks, respectively, per single quark. This can 
be generalized for any baryon with two doubly repre- 
sented quarks D and one singly represented quark S. 
Focusing on the proton and using the charge factors of 
the doubly represented and singly represented quarks as 
2/3 and (-1/3) respectively, we define effective moments 
for single quarks of unit charge, /i^*, by 



1 



Eff 
Ms : 



such that Eq. (|5.7|) becomes 



(5.8) 



(5.9) 



where the charge factors within /i^ 
plicit. 



and are now ex- 
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FIG. 13: Charge radii of A and S* at different quark 
masses. Tfie values for E* are plotted at shifted for clarity. 
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FIG. 15: Charge radii of E*" and H*" at different quark 
masses. The values for E*" are plotted at shifted for 
clarity. 




FIG. 14: Charge radii (magnitude) of E*~ and H*~ at differ- 
ent quark masses. The values for H*~ are plotted at shifted 
for clarity. 



FIG. 16: Ratio of charge radii of the octet/decuplet baryons 
at the ninth quark mass where = 0.215(4) GeV^ 



One could also define a constituent quark mass via 



On the lattice we calculate the baryon magnetic mo- 
ments from the individual quark sector contributions us- 



ing 



11-3) A^5^" 



(5.10) 



In the above equation the factors 2 and 1 in the first 
and second term account for the number of doubly and 
singly represented quarks, while the charges are indicated 
in parentheses. Equating quark sectors in Eqs. (|5.9p and 
()5.10p yields for the effective moments, 



Eff _ q Latt 
— ^•'MS : 

,,Eff _ 3 Tjjtj 



(5.11) 



Eff 

Ms 



Eff _ 

^J■D — 



2m™ 



(5.12) 



revealing that fig^ ~ /i™ in most constituent quark 
models. 

For the decuplet baryons, the magnetic moment is the 
sum of the individual constituent-quark contributions. 
Hence Eq. (|5.9p for the A baryons becomes 



MA+ = 2 ( - ) A^^" + 1 



1 



Eff 

Ms ■ 



(5.13) 



On the lattice this is exactly the equation we use to build 
the decuplet baryon moments from the quark sector con- 
tributions. Therefore, the quark level magnetic moments 
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FIG. 17: Magnetic moments of a u quark in A and H* as 
a function of quark mass. The values for H* are plotted at 
shifted for clarity. 
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FIG. 19: Magnetic moments of an s quark in S* and H* 
for different quark masses. The values for H* are plotted at 
shifted for clarity. 
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FIG. 18: Magnetic moments of a -u quark in E* and H* as 
a function of quark mass. The values for H* are plotted at 
shifted for clarity. 



FIG. 20: Effective moments of the u quark sector in the octet 
and the decuplet baryons at the ninth quark mass where — 
0.215(4) GeV^ 



that we calculate are the effective moments of the quarks 
for both the doubly and singly represented quarks, i.e., 

^Eff ^ ^Latt 

In Figs. [20l and [2T]we plot the effective moments of the 
u and s quark sectors of the octet and decuplet baryons 
at the ninth quark mass. Here we observe that the quarks 
in the octet baryons show far more environmental sensi- 
tivity than their counterparts in the decuplet baryons. 

The baryon magnetic moments are plotted in Figs. 
[23l and [24l For the magnetic moment of fi" we take the 
value of fi^- at the SU(3) limit viz., —1.697± 0.065 fijy, 
which is smaller than the value given by the Particle Data 
Group (-2.02 ±0.05 ^in)- 



This discrepancy could be partly due to the fact that 
the mass of il^ from our lattice calculation (1.73 ± 
0.012 GeV) is slightly larger than the experimentally 
measured value (1.67 GeV). 

Another reason for this discrepancy is likely to reside in 
the absence of KE loops in the virtual decay of il^ . The 
virtual transition Q — > 'E.K requires the presence of a light 
sea-quark flavor, while in QQCD, there is only a heavy 
valence strange quark. In reality this would provide an 
important contribution, since S is a lower mass state than 
O^. The predominant contribution is fJ'^ "EPK^ with 
the z-component of angular momentum in the positive 
direction. This process will act to enhance the magni- 
tude of the negative moment. The absence of such loops 
in QQCD represents missing physics and causes the dis- 
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FIG. 21: Effective moments of the s quark sector in the octet 
and the decuplet baryons at the ninth quark mass where — 
0.215(4) GeV^ 
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FIG. 23: Magnetic moments (magnitude) of E*~ and H*~ 
at different quark masses. The values for H*~ are plotted at 
shifted for clarity. 
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FIG. 22: Magnetic moments of A"'" and E*"^ at different quark 
masses. The values for E*^ are plotted at shifted for 
clarity. 



FIG. 24: Magnetic moments of E*" and H*" at different quark 
masses. The values for H*" are plotted at shifted for 
clarity. 



crepancy from the values of full QCD. This is certainly a 
good place to search for dynamical sea quark effects. 

Figure [25] compares the magnetic moment of the A"*" 
with our earlier result for the proton magnetic moment 
on the same set of gauge field configurations [13] ■ A 
simple quark model predicts that the proton and the 
have equal magnetic moments. Hovifcver the inter- 
play between the different pion-loop contributions to the 
A"*" magnetic moments indicate that the proton magnetic 
moment should be greater than that of the A"*" in Full 
QCD Ml 

The presence of the A Ntt decay channel is par- 
ticularly important for the quark mass dependence of A 
properties [37[. Rapid curvature associated with non- 
analytic behavior is shifted to larger pion masses near 



the A^-A mass splitting, ttItt ~ — Mn- As described 
below, quenched-QCD decay-channel contributions come 
with a sign opposite to that of full QCD. This artifact 
holds tremendous promise for revealing unmistakable sig- 
natures of the quenched meson cloud. 

The change in sign for the decay-channel contributions 
is easily understood through the consideration of the 
quark flow diagrams in Fig. I26[ illustrating the meson- 
cloud contributions to the A++ resonance in full QCD. 
Quark flow diagram (a) corresponds to the hadronic pro- 
cess described in the top left diagram of Fig. [551 Since 
QCD is flavor-blind, the process illustrated in diagram 
(b) is equivalent to diagram (a) provided the masses of 
the u and d quarks are taken to be equal. On its own, 
diagram (b) describes the decay of the A++ to a doubly- 
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FIG. 25: Magnetic moments of and the proton at different 
quark masses. 
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FIG. 26; Quark-fiow diagrams for meson-cloud contributions 
to A++ in full QCD. 



charged uuu "proton," which we denote Of course, 

such states do not exist in full QCD and diagram (c) 
provides a contribution which is exactly equal but oppo- 
site in sign to diagram (b) when the intermediate state 
is a uuu proton. Upon quenching the theory, both dia- 
grams (a) and (b) are eliminated, leaving only diagram 
(c). Hence the physics of the A Ntt decay is present 
in the quenched approximation f33| but its contribution 
has the wrong sign. This signature of quenched chiral 
physics is manifest in our results. 

All the decuplet baryons show, to some extent, a turn 
over of the magnetic moment in the low quark mass re- 
gion just above the opening of the Ntt decay channel. 
The magnitude of the turn over is dampened by the pres- 
ence of a strange quark, which is seen by the fact that 
S* has a smaller turn over than A. The S baryons with 
two s quarks only admit Kaon loops and do not display 
a turn over, further clarifying a link to chiral physics. 

Magnetic moments of the decuplet baryons are listed 
in Tables Em] to KVl 
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TABLE XIII: Magnetic moments of the A"*" in nuclear magnetons for different in GeV^. Quark sector contributions for 
single quarks of unit charge arc also provided. The magnetic moment of A~ is equal in magnitude to that of A"*" with a negative 
sign and that of A*^ is 0. Charge symmetry also requires that the A+^ has a magnetic moment twice that of A+. The magnetic 



moment of Q 


is that of A at the S'?7(3)flavor limit where m% 


= 0.4854(31) GeV^, and takes the value 


-1.697(65) IJ.N- 


ml 


UA (mjv) 


dA (mjv) 


A+ (hm) 


0.9960(56) 


1.466(31) 


1.466(31) 


1.466(31) 


0.8936(56) 


1.512(36) 


1.512(36) 


1.512(36) 


0.7920(55) 


1.559(42) 


1.559(42) 


1.559(42) 


0.6920(54) 


1.604(50) 


1.604(50) 


1.604(50) 


0.6910(35) 


1.615(38) 


1.615(38) 


1.615(38) 


0.5925(33) 


1.662(48) 


1.662(48) 


1.662(48) 


0.4854(31) 


1.697(65) 


1.697(65) 


1.697(65) 


0.3795(31) 


1.749(74) 


1.749(74) 


1.749(74) 


0.2839(33) 


1.792(87) 


1.792(87) 


1.792(87) 


0.2153(35) 


1.80(11) 


1.80(11) 


1.80(11) 


0.1384(43) 


1.79(15) 


1.79(15) 


1.79(15) 


0.0939(44) 


1.60(28) 


1.60(28) 


1.60(28) 



TABLE XIV: 


Magnetic moments of E*, E*'^ 


and E* in 


nuclear magnetons with quark sector contributions for 


a single quark 


of unit charge 


at different values. 










ml 


Me* (/xat) 


se* (mat) 


E*+ (mjv) 




E*- (/xjv) 


0.9960(56) 


1.482(38) 


1.671(51) 


1.418(36) 


-0.063(7) 


-1.545(42) 


0.8936(56) 


1.524(43) 


1.675(54) 


1.474(41) 


-0.050(7) 


-1.574(46) 


0.7920(55) 


1.568(48) 


1.680(57) 


1.531(47) 


-0.037(6) 


-1.605(51) 


0.6920(54) 


1.609(55) 


1.683(62) 


1.585(54) 


-0.025(5) 


-1.634(57) 


0.6910(35) 


1.621(44) 


1.696(48) 


1.596(43) 


-0.025(3) 


-1.646(45) 


0.5925(33) 


1.664(52) 


1.698(54) 


1.653(51) 


-0.011(2) 


-1.675(53) 


0.4854(31) 


1.697(65) 


1.697(65) 


1.697(65) 


0.000(0) 


-1.697(65) 


0.3795(31) 


1.744(72) 


1.709(69) 


1.755(72) 


0.012(2) 


-1.732(71) 


0.2839(33) 


1.781(80) 


1.724(75) 


1.800(82) 


0.019(5) 


-1.762(78) 


0.2153(35) 


1.796(90) 


1.738(80) 


1.816(95) 


0.020(10) 


-1.777(85) 


0.1384(43) 


1.81(11) 


1.766(86) 


1.82(12) 


0.013(21) 


-1.793(97) 


0.0939(44) 


1.76(15) 


1.75(10) 


1.76(18) 


0.004(38) 


-1.75(13) 



TABLE XV: Magnetic 


moments of H*° and H* 


in nuclear magnetons 


with quark sector contributions for a 


single quark of 


unit charge at different 


values. 








ml 


sa* (mjv) 


wh* (mjv) 


H*" (Miv) 


H*- inN) 


0.9960(56) 


1.681(63) 


1.494(48) 


-0.124(18) 


-1.619(57) 


0.8936(56) 


1.683(65) 


1.534(52) 


-0.099(16) 


-1.633(60) 


0.7920(55) 


1.685(67) 


1.576(56) 


-0.073(14) 


-1.649(63) 


0.6920(54) 


1.687(69) 


1.615(62) 


-0.048(11) 


-1.663(67) 


0.6910(35) 


1.697(56) 


1.626(51) 


-0.047(6) 


-1.674(54) 


0.5925(33) 


1.698(59) 


1.666(57) 


-0.021(4) 


-1.687(58) 


0.4854(31) 


1.697(65) 


1.697(65) 


0.000(0) 


-1.697(65) 


0.3795(31) 


1.703(67) 


1.738(69) 


0.023(4) 


-1.714(68) 


0.2839(33) 


1.709(69) 


1.772(73) 


0.042(9) 


-1.730(70) 


0.2153(35) 


1.714(71) 


1.792(78) 


0.052(16) 


-1.740(73) 


0.1384(43) 


1.725(73) 


1.811(86) 


0.057(31) 


-1.753(75) 


0.0939(44) 


1.725(77) 


1.84(10) 


0.076(48) 


-1.763(79) 
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D. Electric Quadrupole Form Factors 

The E2 form factors of the spin-3/2 decuplet baryons 
provide interesting information about the distribution of 
charge and its deviation from spherical symmetry. In 
Fig. [27] we show the correlation function proportional to 
the E2 quadrupole form factor of the u quark in A at 
the 5'?7(3)fiavor limit, in units of e/M^, as a function of 
Euclidean time. Figure indicates the quadrupole form 
factor of the u quark in A at the ninth quark mass. Here 
the employment of the splittings technique facilitates the 
extraction of the signal. In both cases a nontrivial result 
is obtained. As mentioned in Sec. Ill Dl we consider the 
symmetry of the last two terms in Eq. (j2.35p as the decid- 
ing factor in selecting the upper limit of the fit-window. 

The quark sector contributions to the form factors in 
units of e/Mjj of all the decuplet members are indicated 
in Tables IXVIl to IXVIIIl For an axially deformed ob- 
ject the quadrupole form factor is related to the charge 
distribution in the Brcit frame through Q 



0.4 



fe2(0) = Ml / dVV;(r)(3z2 _ r2)^(r) , 



(5.15) 



where 3z^ — is the standard operator used for 
quadrupole moments. A positive quadrupole form factor 
for a positively charged baryon indicates a prolate charge 
distribution, while a negative quadrupole form factor in- 
dicates an oblate charge distribution. In non-relativistic 
models, the E2 form factor vanishes unless some configu- 
ration mixing of higher orbital-angular momentum states 
is included in the baryon ground state. 

The E2 form factors of the charged decuplet baryons 
in units of fm^ for different values of are listed in 
table IXiXl The E2 form factor of the A° is identically 
equal to zero and for the other neutral baryons, it is close 
to zero. The results for the charged decuplet baryons are 
non-zero, indicating that they have a deformed shape. 

The quark sector contributions to the E2 form factors 
are shown in Figs. [25] to [HH Once again, the importance 
of chiral physics is manifest in these results. A significant 
enhancement of the magnitude of the light-quark sector 
contribution to the E2 form factor is observed in the A 
as the opening of the Nn decay channel is approached. A 
similar effect is seen, to a lesser extent, in the E*, while 
for the 'Er no chiral curvature is observed due to the two 
s quarks admitting only Kaon loops, as discussed in the 
previous section. The E2 form factors for the various 
decuplet baryons are shown in Figs. [5^ to [Ml 

From our simulation we conclude that the E2 form fac- 
tor of the 0~ baryon (the value of the A^ form factor 
at 5t/(3)fiavor limit) is (0.86 ±0.12) x 10^^ fj^2_ rpj^g 
curacy of our result indicates a definite non-zero value of 
the E2 form factor of il^, and we favor a positive value. 
Since f2^ is a negatively charged baryon, this result im- 
plies that it has an oblate shape, with the equatorial axis 
being larger than the polar axis. 

Similarly the E2 form factor for the A"*" is (—0.86 ± 
0.12) X 10~^ fm^ at the 51/(3) flavor symmetry point. 




16 18 
Euclidean Time 



FIG. 27: E2 electric form factor of the u or d quark sector of 
the A at the 5'il7(3)flavor limit as a function of time. 




16 18 20 

Euclidean Time 
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FIG. 28: E2 electric form factor of the m or d quark sector 
of the A as a function of time at the ninth quark mass where 
ml = 0.215(4)GeV^ 



Our results for the A+ compare favorably with the re- 
sults of Ref . 8] . Using the closest available pion masses 
of 533(3) MeV, and 563(4) MeV for this study, and [| re- 
spectively, we find the E2 form factor to be (— 1.16±29) x 
10-2 fm^, to be compared with (-1.08 ±40) x 10"^ fm^ 
We note however, that this study is performed at a finite 
Q2 ^ 0.230(1) GeV^ whilst [1] reports results at = q. 

The negative E2 form factor of a positive A+ baryon 
implies an oblate shape in the Breit frame. As illustrated 
in Fig. [321 the E2 form factor grows substantially in mag- 
nitude as the chiral limit is approached, taking the value 
—0.030(11) fm^ at our lightest quark mass. We note that 
the E2 form factor of A++ is twice that of the A+ E2 
form factor and hence takes the value 
our lightest quark mass. 



-0.060(23) fm2 at 
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TABLE XVI: Quark sector contributions to the E2 form factor of A at = 0.230(1) GeV^ in fixed units of e/m%. Sector 
contributions are for a single quark having unit charge. The fit windows are selected using the criteria outlined in Ref. 14]. 



mi (GeV^) 
















fit V3,lU6 


fit window 


Adof 


fit V3,lU6 


fit window 


Adof 


0.9960(56) 


—0.117(16) 


16 — 20 


1.63 


—0.117(16) 


16 — 20 


1.63 


0.8936(56) 


—0.123(18) 


16 — 20 


1.39 


—0.123(18) 


16 — 20 


1.39 


0.7920(55) 


—0.130(22) 


16 — 20 


1.16 


—0.130(22) 


16 — 20 


1.16 


0.6920(54) 


-0.137(26) 


16 - 20 


1.01 


-0.137(26) 


16 - 20 


1.01 


0.6910(35) 


-0.163(17) 


16 — 20 


0.76 


-0.163(17) 


16 — 20 


0.76 


0.5925(33) 


-0.177(21) 


16 - 20 


0.75 


-0.177(21) 


16 - 20 


0.75 


0.4854(31) 


-0.194(27) 


16 - 20 


0.86 


-0.194(27) 


16-20 


0.86 


0.3795(31) 


-0.218(40) 


16-19 


1.03 


-0.218(40) 


16-19 


1.03 


0.2839(33) 


-0.263(67) 


16-19 


1.57 


-0.263(67) 


16-19 


1.57 


0.2153(35) 


-0.32(11) 


16 - 19 


1.20 


-0.32(11) 


16-19 


1.20 


0.1384(43) 


-0.52(20) 


16-18 


0.72 


-0.52(20) 


16-18 


0.72 


0.0939(44) 


-0.68(26) 


15 - 16 


1.06 


-0.68(26) 


15 - 16 


1.06 



TABLE XVII: Quark sector contributions to the E2 form factor of E* at = 0.230(1) GeV^ in units of e/m%. Sector 
contributions are for a single quark having unit charge. The fit windows are selected using the criteria outlined in Ref. [l3 |. 



ml (GeV^) 




Its* or ds* 






St.* 






fit value 


fit window 


Xdof 


fit value 


fit window 


Xdof 


0.9960(56) 


-0.132(19) 


16 - 20 


1.52 


-0.113(29) 


16-20 


0.75 


0.8936(56) 


-0.136(22) 


16 - 20 


1.36 


-0.118(31) 


16-20 


0.71 


0.7920(55) 


-0.145(25) 


16-20 


1.20 


-0.124(33) 


16-20 


0.66 


0.6920(54) 


-0.147(29) 


16 - 20 


1.07 


-0.131(36) 


16-20 


0.65 


0.6910(35) 


-0.172(18) 


16-20 


0.91 


-0.170(22) 


16-20 


0.48 


0.5925(33) 


-0.183(22) 


16-20 


0.89 


-0.180(24) 


16-20 


0.61 


0.4854(31) 


-0.194(27) 


16 - 20 


0.86 


-0.194(27) 


16-20 


0.86 


0.3795(31) 


-0.208(36) 


16-20 


0.48 


-0.211(32) 


16-20 


1.89 


0.2839(33) 


-0.225(51) 


16-17 


0.41 


-0.231(38) 


16-17 


0.61 


0.2153(35) 


-0.233(73) 


16 - 19 


1.08 


-0.257(48) 


16 - 19 


1.84 


0.1384(43) 


-0.29(11) 


16-17 


1.71 


-0.300(67) 


16 - 17 


1.07 


0.0939(44) 


-0.42(16) 


16-17 


0.94 


-0.325(88) 


16-17 


0.31 



E. Magnetic Octupole Moments 

The magnetic octupole form factors are calculated on 
the lattice by considering a combination of ratios of three 
and two-point functions as given in Eq. (|2.32p . 

Figure [551 provides a plot of the correlator proportional 
to the M3 form factor of a u quark in the A as a function 
of Euclidean time at the SU{?>) flavor limit. Figure [36] 
provides the M3 form factor at the ninth quark mass, 
where a plateau is realized using the splittings method. 
Tables [XXl to IXXTIl list the quark sector Mi form factors. 
The magnetic octupole form factors (M3) of the decuplet 
baryons are listed in Table [XXIIII 

Plots of the quark sector contributions to the Af3 form 
factors are provided in Figs. [37] to [351 Figs. [301 to [151 show 



the M3 form factors for the decuplet baryons. 

Like E2, the Mi form factors require nonzero orbital 
angular momentum admixtures in the ground state wave 
function [?[. Our statistics are sufficient to reveal a non- 
trivial result for the M3 form factor of the A for the 
first time. We find a result of —3.7(3.3) e/2m'^ for the 
A+ at a squared pion mass of 0.094(4) GeV^, close to the 
physical limit. At larger masses the results are consistent 
with zero, but systematically negative, with enhancement 
at the lightest masses for the A, and to a lesser extent 
the E*. 
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TABLE XVIII: Quark sector contributions to the E2 form factor of H* at = 0.230(1) GeV^ in units of e/m%. Sector 
contributions are for a single quark having unit charge. The fit windows are selected using the criteria outlined in Ref. 14]. 



(GeV^) 




Ss* 






Us* or (is* 






nt value 


lit window 


Xdof 


nt value 


nt window 


,,2 
Xdof 


0.9960(56) 


-0.127(36) 


16 - 20 


0.64 


-0.157(25) 


16-20 


1.75 


0.8936(56) 


-0.131(37) 


16 - 20 


0.64 


-0.157(27) 


16-20 


1.52 


0.7920(55) 


-0.136(39) 


16 - 20 


0.64 


-0.159(30) 


16-20 


1.35 


0.6920(54) 


-0.141(40) 


16-20 


0.66 


-0.159(33) 


16-20 


1.18 


0.6910(35) 


-0.180(24) 


16-20 


0.65 


-0.184(20) 


16-20 


1.17 


0.5925(33) 


-0.186(25) 


16-20 


0.73 


-0.190(23) 


16-20 


1.04 


0.4854(31) 


-0.194(27) 


16-20 


0.86 


-0.194(27) 


16-20 


0.86 


0.3795(31) 


-0.201(29) 


16-17 


0.82 


-0.198(33) 


16 - 21 


0.62 


0.2839(33) 


-0.208(31) 


16-17 


0.59 


-0.200(41) 


16-17 


0.23 


0.2153(35) 


-0.214(34) 


16-17 


0.80 


-0.191(52) 


16-17 


0.06 


0.1384(43) 


-0.222(38) 


16-17 


0.92 


-0.184(71) 


16-18 


0.34 


0.0939(44) 


-0.222(41) 


15-16 


0.17 


-0.183(82) 


15-16 


0.70 



TABLE XIX: E2 form factors at = 0.230(1) GeV^ of the charged decuplet baryons in units of lO^^fm^ for different ml 
values. The E2 form factor of the A~ at the S'[/(3)flavor limit where = 0.485(3) provides the E2 form factor of Q~ . 



ml (GeV^) 


A++ 


A+ 


A- 








0.9972(55) 


-1.03(14) 


-0.517(69) 


0.517(69) 


-0.613(77) 


0.555(97) 


0.60(14) 


0.8936(56) 


-1.08(16) 


-0.541(80) 


0.541(80) 


-0.629(88) 


0.57(11) 


0.62(15) 


0.7920(55) 


-1.15(19) 


-0.575(96) 


0.575(96) 


-0.65(10) 


0.60(12) 


0.63(16) 


0.6920(54) 


-1.21(23) 


-0.61(12) 


0.61(12) 


-0.67(12) 


0.63(14) 


0.65(17) 


0.6910(35) 


-1.44(15) 


-0.718(75) 


0.718(75) 


-0.765(77) 


0.757(86) 


0.80(10) 


0.5925(33) 


-1.56(18) 


-0.782(91) 


0.782(91) 


-0.813(92) 


0.804(99) 


0.83(11) 


0.4854(31) 


-1.71(24) 


-0.86(12) 


0.86(12) 


-0.86(12) 


0.86(12) 


0.86(12) 


0.3795(31) 


-1.93(35) 


-0.96(18) 


0.96(18) 


-0.91(17) 


0.92(15) 


0.88(13) 


0.2839(33) 


-2.32(59) 


-1.16(29) 


1.16(29) 


-0.99(25) 


1.00(20) 


0.91(15) 


0.2153(35) 


-2.79(95) 


-1.40(48) 


1.40(48) 


-0.99(37) 


1.07(28) 


0.91(17) 


0.1384(43) 


-4.6(1.8) 


-2.31(88) 


2.31(88) 


-1.28(58) 


1.30(41) 


0.92(21) 


0.0939(44) 


-6.0(2.3) 


-3.0(1.1) 


3.0(1.1) 


-1.99(86) 


1.71(57) 


0.92(23) 
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TABLE XX: Quark sector contributions to the M3 form factor of A at = 0.230(1) GeV^ in units of e/2m%. Sector 
contributions are for a single quark having unit charge. The fit windows are selected using the criteria outlined in Rcf. 14]. 



(GeV^) 
















fit value 


fit window 


Ador 


fit value 


fit window 


Adoi 


0.9972(55) 


0.19(58) 


18 — 22 


0.83 


0.19(58) 


18 — 22 


0.83 


0.8936(56) 


0.19(68) 


18 — 22 


0.63 


0.19(68) 


18 — 22 


0.63 


0.7920(55) 


0.18(82) 


18 - 22 


0.50 


0.18(82) 


18 - 22 


0.50 


0.6920(54) 


0.2(1.0) 


18 - 22 


0.44 


0.2(1.0) 


18 - 22 


0.44 


0.6910(35) 


-0.10(75) 


18 - 22 


0.74 


-0.10(75) 


18 - 22 


0.74 


0.5925(33) 


-0.09(87) 


18-20 


1.01 


-0.09(87) 


18-20 


1.01 


0.4854(31) 


-0.2(1.2) 


18-20 


0.92 


-0.2(1.2) 


18-20 


0.92 


0.3795(31) 


-0.4(1.4) 


16- 18 


0.44 


-0.4(1.4) 


16- 18 


0.44 


0.2839(33) 


-0.8(1.6) 


15 - 17 


0.01 


-0.8(1.6) 


15 - 17 


0.01 


0.2153(35) 


-1.2(1.9) 


15 - 17 


0.12 


-1.2(1.9) 


15 - 17 


0.12 


0.1384(43) 


-2.4(2.5) 


15 - 17 


0.37 


-2.4(2.5) 


15 - 17 


0.37 


0.0939(44) 


-3.7(3.3) 


15 - 17 


0.32 


-3.7(3.3) 


15 - 17 


0.32 



TABLE XXI: Quark sector contributions to the M3 form factor of E* baryons at — 0.230(1) GeV^ in units of e/2m%. Sector 
contributions are for a single quark having unit charge. The fit windows are selected using the criteria outlined in Ref. fT3]. 



ml (GeV^) Its* or ds* se* 





fit value 


fit window 


Xdof 


fit value 


fit window 


Xdof 


0.9972(55) 


0.13(81) 


18- 


- 22 


0.56 


0.64(80) 


18 - 


- 22 


0.50 


0.8936(56) 


0.12(91) 


18- 


- 22 


0.50 


0.63(90) 


18 - 


- 22 


0.44 


0.7920(55) 


0.1(1.0) 


18- 


- 22 


0.48 


0.6(1.0) 


18 - 


- 22 


0.39 


0.6920(54) 


0.1(1.2) 


18- 


- 22 


0.47 


0.6(1.2) 


18 - 


- 22 


0.38 


0.6910(35) 


-0.16(89) 


18- 


- 22 


0.96 


0.10(89) 


18 - 


- 22 


0.64 


0.5925(33) 


-0.1(1.0) 


18- 


-20 


1.03 


-0.1(1.0) 


18 - 


- 20 


0.90 


0.4854(31) 


-0.2(1.2) 


18- 


- 20 


0.92 


-0.2(1.2) 


18 - 


- 20 


0.92 


0.3795(31) 


-0.3(1.3) 


16- 


- 18 


0.43 


-0.3(1.3) 


16- 


- 18 


0.61 


0.2839(33) 


-0.5(1.4) 


15 - 


- 17 


0.01 


-0.4(1.4) 


15 - 


- 17 


0.30 


0.2153(35) 


-0.7(1.6) 


15 - 


- 17 


0.01 


-0.5(1.5) 


15- 


- 17 


0.38 


0.1384(43) 


-1.1(1.8) 


15 - 


- 17 


0.21 


-0.8(1.7) 


15- 


- 17 


0.07 


0.0939(44) 


-1.2(2.0) 


15 - 


- 17 


0.07 


-1.1(2.0) 


15 - 


- 17 


0.89 
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TABLE XXII: Quark sector contributions to the M3 form factor of H* baryons at = 0.230(1) GeV^ in units of e/2m%. 
Sector contributions are for a single quark having unit charge. The fit windows are selected using the criteria outlined in 
Ref. [il. 



ml (GeV^) SH* Its* or ds* 





fit value 


fit window 


Xdof 


fit value 


fit window 


Xdof 


0.9972(55) 


0.6(1.1) 


18 


- 22 


0.39 


-0.1(1.2) 


18- 


22 


0.45 


0.8936(56) 


0.6(1.2) 


18 


- 22 


0.39 


-0.1(1.2) 


18- 


22 


0.48 


0.7920(55) 


0.5(1.3) 


18 


- 22 


0.40 


-0.1(1.3) 


18- 


22 


0.52 


0.6920(54) 


0.5(1.4) 


18 


- 22 


0.41 


0.0(1.4) 


18- 


22 


0.53 


0.6910(35) 


-0.2(1.1) 


18 


- 22 


0.83 


-0.3(1.1) 


18- 


22 


1.27 


0.5925(33) 


-0.2(1.1) 


18 


-20 


0.92 


-0.2(1.1) 


18- 


20 


1.03 


0.4854(31) 


-0.2(1.2) 


18 


- 20 


0.92 


-0.2(1.2) 


18- 


20 


0.92 


0.3795(31) 


-0.3(1.2) 


16 


- 18 


0.62 


-0.3(1.2) 


16- 


18 


0.54 


0.2839(33) 


-0.3(1.3) 


15 


- 17 


0.52 


-0.4(1.3) 


15 - 


17 


0.17 


0.2153(35) 


-0.3(1.3) 


15 


- 17 


0.63 


-0.4(1.4) 


15 - 


17 


0.05 


0.1384(43) 


-0.3(1.4) 


15 


- 17 


0.23 


-0.6(1.4) 


15 - 


17 


0.47 


0.0939(44) 


-0.4(1.4) 


15 


- 17 


0.44 


-0.3(1.5) 


15 - 


17 


0.82 



TABLE XXIII: M3 form factor results at = 0.230(1) GeV^ of the charged decuplet baryons in units of e/2m% for different 
values. The M3 form factor of the A~ at the S'(7(3)flavor limit where — 0.485(3) provides the M3 form factor of QT . 



ml (GeV^) 


A++ 


A+ 


A" 








0.9972(55) 


0.4(1.2) 


0.19(58) 


-0.19(58) 


-0.04(86) 


-0.30(79) 


-0.4(1.1) 


0.8936(56) 


0.4(1.4) 


0.19(68) 


-0.19(68) 


-0.05(94) 


-0.29(89) 


-0.4(1.2) 


0.7920(55) 


0.4(1.6) 


0.18(82) 


-0.18(82) 


-0.1(1.1) 


-0.3(1.0) 


-0.3(1.3) 


0.6920(54) 


0.4(2.0) 


0.2(1.0) 


-0.2(1.0) 


-0.1(1.2) 


-0.3(1.2) 


-0.3(1.4) 


0.6910(35) 


-0.2(1.5) 


-0.10(75) 


0.10(75) 


-0.18(90) 


0.14(88) 


0.2(1.1) 


0.5925(33) 


-0.2(1.7) 


-0.09(87) 


0.09(87) 


-0.2(1.0) 


0.1(1.0) 


0.2(1.1) 


0.4854(31) 


-0.5(2.4) 


-0.2(1.2) 


0.2(1.2) 


-0.2(1.2) 


0.2(1.2) 


0.2(1.2) 


0.3795(31) 


-0.8(2.8) 


-0.4(1.4) 


0.4(1.4) 


-0.3(1.3) 


0.3(1.3) 


0.3(1.2) 


0.2839(33) 


-1.6(3.1) 


-0.8(1.6) 


0.8(1.6) 


-0.6(1.4) 


0.5(1.4) 


0.3(1.3) 


0.2153(35) 


-2.4(3.8) 


-1.2(1.9) 


1.2(1.9) 


-0.7(1.6) 


0.6(1.5) 


0.3(1.3) 


0.1384(43) 


-4.9(4.9) 


-2.4(2.5) 


2.4(2.5) 


-1.2(1.9) 


1.0(1.7) 


0.4(1.4) 


0.0939(44) 


-7.4(6.6) 


-3.7(3.3) 


3.7(3.3) 


-1.3(2.2) 


1.2(1.9) 


0.4(1.4) 



TABLE XXIV: Collected results for the SI . Results are obtained from the A at the S'C/(3)flavor limit, where — 
0.4854(31) GeV^ 



quantity fit value 



mass (GeV) 1.732(12) 

charge radius (fm^) -0.307(15) 

magnetic moment (/ijv) —1.697(65) 

E2 form factor (10"^fm^) 0.86(12) 

Af3 form factor {e/2m%) 0.2(1.2) 
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FIG. 29: E2 form factor contributions from the u quark sec- FIG. 31: E2 form factor contributions from the s quark sec- 
tors of the A and H*. tors of the S* and H*. 
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FIG. 30: E2 form factor contributions from the u quark sec- 
tors of the E* and H* . 



FIG. 32: Values of the E2 form factors in units of lO'^fm^ 
for the and E*'*' at different quark masses. The values for 
A+ are plotted at shifted for clarity. 
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FIG. 35: M3 form factor of the A at the S'[/(3)fla 
a function of time. 
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FIG. 34: E2 form factors in units of 10"^fm^ for the E*° and 
H*° at different quark masses. The values for are plotted 
at shifted for clarity. 



FIG. 36: A/3 form factor (splitting) of the u quark sector of 
the A at the ninth quark mass where — 0.215(4)GeV^ as 
a function of time. 



28 



T r 



T 1 1 r 



1 r 



J I I I I I I I I L 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 
m ^ (GeV^) 



5 

C\2 



o 

o 
to 



o 

CO 



-2 



-8 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 
m ^ (GeV^) 



FIG. 37: A/3 form factor contributions from tiie u quark sec- 
tors of the A and E*. The results for the E* are offset for 
clarity. 



FIG. 39: M3 form factor contributions from the s quark sec- 
tors of the E* and H*. The results for the H* have been offset 
for clarity. 
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FIG. 38: M3 form factor contributions from the u quark sec- 
tors of the E* and H*. The resuhs for the H* have been offset 
for clarity. 



FIG. 40: M3 form factors for the A+ and S*+ at different 
quark masses. The values for E*+ are plotted at shifted 
for clarity. 



VI. SUMMARY 

We have performed an extensive calculation of the 
electromagnetic properties of decuplet baryons at both 
the quark level and the baryon level, including the 
quadrupole and octupole form factors of the spin-3/2 
baryons. For the first time we obtain non-trivial results 
for both the E2 and M3 form factors. In particular, we 
find decuplet baryons to be oblate in shape. 

We find that the quarks in the decuplet are not as 
sensitive to their environment as their octet counterparts. 
Of particular note, is the discovery that the decuplet- 
baryon radii are smaller than that of the octet baryons. 



contradicting the simple quark model, but substantiating 
hints in the early study of Rcf. 7] . 

A particularly interesting finding is that the suppres- 
sion of sea-quark loop contributions in QQCD reduces 
the decuplet magnetic moment considerably, resulting in 
a turnover in the magnetic moment at light quark masses, 
as illustrated in Fig. [551 At large pion masses, the A+ 
moment is enhanced relative to the proton moment in 
accord with earlier quenched lattice QCD calculations 
0, and model expectations. However, as the chiral 
regime is approached, the non-analytic behavior of the 
quenched meson cloud is revealed, enhancing the proton 
and suppressing the , in accord with the expectations 
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FIG. 41: M3 form factors for the E*~ and H*~ at different 
quark masses. The values for H*~ are plotted at shifted 
for clarity. 



of QxPT. This suppression should be absent in full QCD. 
Hence this is one particular case that we have identified 
as a place to look for effects of unquenching in future dy- 
namical simulations. We also predict that unquenching 
effects should be observed in the magnetic moment, 
which we find to be suppressed in quenched QCD with 
regard to the experimentally measured value due to the 
absence of KE loops in the virtual decay of . 

Through a calculation of the decuplet E2 form fac- 
tors and electric quadrupole moments, we predict oblate 
shapes for the decuplet baryons. It will be interesting to 
confront this prediction with an experimental measure- 
ment of the Q~ quadrupole form factor. We provide a 
summary of all the fl^ measurements in Table IXXIVI 

Finally we have obtained non-trivial values for the Af 3 
form factor for the first time in lattice QCD studies. 
These results provide an interesting and novel forum for 
the further development of our understanding of non- 
perturbative QCD. 



1.0 



0.5 



0.0 



u 
o 

o -0.5 
cd 

B -1.0 

o 

„ -1.5 
-2.0 



T 1 1 1 1 1 1 1 r 



1| i\ ^ " i 



J I I I I I I I I L 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 
m = (GeV^) 
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